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TITLE OF THE INVENTION 


AUTOMATED METHOD AND SYSTEM FOR THE SEGMENTATION OF LUNG 
REGIONS IN COMPUTED TOMOGRAPHY SCANS 

CROSS-REFERENCE TO PROVISIONAL APPLICATION 
This application claims the benefit of and priority to U.S. Provisional Application 
Serial No. 60/176,297, filed January 18, 2000. 

STATEMENT REGARDING FEDERALLY SPONSORED RESEARCH 

The U.S. Government has a paid-up license in this invention and the rights in limited 
circumstances to require the patent owner to license others on reasonable terms as provided 
for by the terms of grant numbers CA48985, CA62625, and CA64370 awarded by USPHS. 

BACKGROUND OF THE INVENTION 

Field of the Invention: 

The invention relates generally to a method and system for the computerized, fully- 
automated delineation of the lung regions in thoracic computed tomography (CT) images. 
Novel developments and implementations include techniques for elimination of the patient 
table from the segmented thoracic region, segmentation of the trachea and main bronchi to 
prevent their inclusion within the segmented lung regions, a two- and three-dimensional 
rolling ball algorithm to further refine the segmented lung regions, separation of the two 
lungs at the anterior junction line, and identification of the diaphragm to prevent its inclusion 
within the segmented lung regions. 

The present invention also generally relates to computerized techniques for automated 
analysis of digital images, for example, as disclosed in one or more of U.S. Patents 
4,839,807; 4,841,555; 4,851,984; 4,875,165; 4,907,156; 4,918,534; 5,072,384; 5,133,020; 
5,150,292; 5,224,177; 5,289,374; 5,319,549; 5,343,390; 5,359,513; 5,452,367; 5,463,548; 
5,491,627; 5,537,485; 5,598,481; 5,622,171; 5,638,458; 5,657,362; 5,666,434; 5,673,332; 
5,668,888; 5,740,268; 5,790,690; 5,832,103; 5,873,824; 5,881,124; 5,931,780; 5,974,165; 
5,982,915; 5,984,870; 5,987,345; and 6,011,862; as well as U.S. patent applications 


08/173,935; 08/398,307 (PCT Publication WO 96/27846); 08/536,149; 08/562,087; 
08/900,188; 08/900,189; 08/900,191; 08/900,361; 08/979,623; 08/979,639; 08/982,282; 
09/027,468; 09/027,685; 09/028,518; 09/053,798; 09/092,004; 09/121,719; 09/131,162; 
09/141,535; 09/156,413; 09/298,852; and 09/471,088; PCT patent applications 
5 PCT7US99/24007; PCT/US99/25998; and U.S. provisional patent applications 60/160,790 
and 60/176,304; all of which are incorporated herein by reference. 

The present invention includes the use of various technologies referenced and 
described in the above-noted U.S. Patents and Applications, as well as described in the 
references identified in the appended APPENDIX and cross-referenced throughout the 
10 specification by reference to the corresponding number, in brackets, of the respective 
I) references listed in the APPENDIX, the entire contents of which, including the related 

patents and applications listed above and the references listed in the APPENDIX, are 
-\ incorporated herein by reference. 

? 1 5 Discussion of the Background: 

Helical computed tomography (CT) of the thorax is widely used to evaluate numerous 
lung diseases, including lung nodules, emphysema, and pulmonary embolism [1]. The recent 

^ availability of multi-slice CT scanners promises to expand further the role of CT in the 

diagnostic imaging arena. The increasing volume of thoracic CT studies and the concomitant 

20 burgeoning of image data these studies generate have prompted many investigators to 

develop computer-aided diagnostic (CAD) methods to assist radiologists in evaluating CT 
images [2-12]. To provide radiologists with useful and reliable information, most such CAD 
methods will require accurate identification of the lung boundaries within the images, a 
preprocessing step known as "lung segmentation." 

25 The requirement for accurate lung segmentation is two-fold. First, the pathologies 

that continue to motivate development of the aforementioned CAD schemes are 
predominantly located within or impact the lungs. Consequently, these CAD schemes are 
designed to accommodate the anticipated appearance of the lung regions in CT images. 
Moreover, spatially limiting further processing to the lungs greatly reduces computation time, 

30 since the lungs occupy a fraction of the total volume of data acquired during a CT scan. 

Second, lung segmentation must be complete since abnormalities such as lung nodules may 
exist at the extreme periphery of the lungs. If the entire lung is not segmented, such 
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abnormalities will be lost to subsequent analyses. Moreover, quantitative assessment of lung 
volume for the evaluation of, for example, emphysema, will be compromised by erroneous 
lung segmentation. 

Aside for its application as a preprocessing step for CAD methods, automated lung 
5 segmentation may be useful for image data visualization. The three-dimensional display of 
CT image data is an area of rapid development with a number of well documented clinical 
applications [13]. Initial lung segmentation would be required in a situation where, for 
example, a volume-rendered version of the lung parenchyma is desired as an aid to the 
radiologist's diagnostic task. 
1 0 Some form of automated lung segmentation serves as a necessary preprocessing step 

for various CAD schemes [2,3,5,14-26]. However, there is a need to improve and refine 
known techniques for automated lung segmentation. 

" SUMMARY OF THE INVENTION 

15 Accordingly, an object of this invention is to provide an improved method and system 

for segmenting the lung regions in thoracic CT images. 

It is another object of this invention to provide an automated method and system for 
segmenting the thorax region and eliminating the patient table from thoracic CT images. 

It is a further object of this invention to provide an automated method and system for 
20 delineating the trachea and main bronchi within the segmented thorax to prevent these 
structures from erroneously contributing to the segmented lung regions. 

It is yet another object of this invention to provide an automated method and system 
for delineating the anterior junction line to prevent the merging of left and right lungs as a 
single segmented lung region. 
25 It is a still further object of this invention to provide an automated method and system 

for refining the segmented lung regions through a two- and three-dimensional rolling ball 
algorithm to ensure proper inclusion of pixels along the peripheral and mediastinal aspects of 
the lungs. 

It is still yet a further object of this invention to provide an automated method and 
30 system for identifying the diaphragm region within the thorax to prevent the erroneous 
inclusion of the diaphragm within the segmented lung regions. 
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These and other objects are achieved according to the invention by providing a new 
and improved automated method, storage medium storing a program for performing the steps 
of the method, and system in which segmentation of lung regions within thoracic images is 
performed. The method, on which the system is based, includes the steps of acquiring image 
5 data representative of a cross-sectional thoracic image; generating initial lung contours to 
segment the lung regions in the cross-sectional thoracic image; identifying within the lung 
region at least one portion corresponding to the diaphragm; and excluding from the lung 
regions the at least one portion corresponding to the diaphragm. 

Preferably the step of identifying includes identifying holes within the lung regions; 
10 determining, for each hole, a geometric feature; comparing the geometric feature of each hole 
* with a threshold; and determining, for each hole, whether the hole corresponds to the 

I diaphragm based on the comparison of the geometric feature to the threshold; and the step of 

excluding step includes excluding from the lung regions the holes corresponding to the 
diaphragm. 

1 5 The method also includes identifying an anterior junction line; extracting from the 

lung regions pixels along the anterior junction line to separate the lung regions; identifying 
within the lung regions portions corresponding to the trachea and main stem bronchi; 
excluding from the lung regions the portions corresponding to the trachea and the main stem 
bronchi; refining the lung contours by applying a rolling ball filter to the lung contours to 

20 identify indentations along the lung contours; determining, for each indentation identified by 
the rolling ball filter, whether the indentation corresponds to the diaphragm; and preventing 
the rolling ball filter from including within the segmented lung regions the indentations 
corresponding to the diaphragm. Preferably the rolling ball filter is a three-dimensional 
rolling ball filter applied to the lung contours in the cross-sectional thoracic image and to 

25 other lung contours in other cross-sectional thoracic images. 

Specific application is given for the delineation of the lung regions in standard helical 
CT scans for such computerized detection schemes such as lung nodule detection, 
emphysema detection, and pulmonary embolism detection. For example, all or portions of 
the segmentation scheme of the present invention may be applied in the three-dimensional 

30 rendering of lung regions and/or in the detection of lung nodules as described in U.S. 

provisional patent application serial number 60/176,304. However, the described methods 
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are also valid for images acquired by conventional CT, high resolution CT, multi-slice CT, 
and low-dose helical CT. 

BRIEF DESCRTPTTON OF THE DRAWINGS 
5 A more complete appreciation of the invention and many of the attendant advantages 

thereof will be readily obtained as the same becomes better understood by reference to the 
following detailed description when considered in connection with the accompanying 
drawings, wherein: 

Figure 1 is a flowchart illustrating a method for the automated segmentation of the 
10 lung regions in thoracic CT images; 

Figure 2 A is an image illustrating the selection of a gray level threshold for thorax 
Ij: segmentation based on a cumulative gray level histogram obtained from pixels along a 
i, diagonal in the image; 

Figure 2B is the cumulative gray level histogram obtained from the pixels along the 
1 5 diagonal line in the image of Figure 2A; 

I Figures 3 A and 3B depict the appearance of the segmented thorax region prior to, and 

subsequent to, respectively, the elimination of portions of the patient table that might be 
"connected" to the thorax region after initial segmentation; 

Figures 4A and 4B are flowcharts illustrating a method for segmentation of the 
20 trachea and the main stem bronchi; 

Figures 5 A and 5B are images demonstrating the results of trachea and main bronchi 
segmentation, respectively, within the segmented thorax regions; 

Figure 6A is a typical gray level histogram constructed from pixels within the 
segmented thorax region; 
25 Figure 6B is an image illustrating the results of initial lung segmentation; 

Figure 7 is a flowchart illustrating a method for separating the right and left lungs 
merged at the anterior junction line; 

Figure 8 A is an image of the anterior junction line and an identified cleft point; 
Figure 8B is the image of Figure 8 A with the pixels eliminated along the delineated 
30 anterior junction line; 

Figure 9 is an image depicting the exclusion of dense structures such as juxta-pleural 
nodules and hilar vessels from the initial lung segmentation contours; 
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Figures 10A and 10B are schematic illustrations demonstrating the two-dimensional 
rolling ball algorithm applied to the external aspect of the initial lung segmentation contours 
to identify and appropriately rectify erroneous indentations in the contours; 

Figures 1 1 A and 1 IB are images of segmented lung regions before and after, 
5 respectively, application of the rolling ball algorithm; 

Figures 12 A is a lung region binary image with a large, circular hole (identified with 
an arrow) caused by the diaphragm; 

Figure 12B is a lung region image demonstrating how the rolling ball algorithm may 
incorrectly include an indentation caused by the diaphragm within the segmented lung 
10 region; 

I Figures 13 A and 13B are flowcharts illustrating two methods for identifying pixels 

» , that belong to the diaphragm and excluding such pixels from the segmented lung regions; 
=. Figure 14 is a block diagram illustrating a system for implementing the inventive 

method for segmenting lung regions in thoracic CT images; and 
"1 5 Figure 15 is a schematic illustration of a general purpose computer system 1500 

- : programmed according to the teachings of the present application. 

I DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

Referring now to the drawings, and more particularly to Figure 1 thereof, a flowchart 

20 of the automated method for the segmentation of the lung regions in thoracic CT scans is 

shown. The overall scheme includes an initial acquisition of CT image data. A database of 
17 helical thoracic CT cases (493 individual section images) was used to develop this method 
and evaluate the performance of various algorithms. The CT section images in the database 
were 5 12 by 5 12 pixels. 

25 In step 101 the thorax is segmented from the background image data. Specifically, a 

cumulative gray level profile is constructed from pixels along the diagonal of each CT 
section image, and the shape of this profile is used to identify a gray level threshold. The 
gray level threshold is applied to the image so that the brightest pixels in the image remain 
"on" in the binary image that is created. A contour-detection algorithm is used to identify the 

30 outer margin of the largest "on" region in the binary image, and the set of all image pixels 
that lie within this contour is considered the segmented thorax region. 
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In step 103 trachea and main bronchi are segmented in all sections in which they 
appear. Pixels identified as belonging to the trachea or main bronchi are effectively 
eliminated from the segmented thorax region to prevent subsequent inclusion within the 
segmented lung regions. 
5 Initial lung segmentation is performed in step 105. A gray level histogram is 

constructed from the remaining pixels in the segmented thorax region, and the broad 
minimum between the peaks of this typically bimodal histogram is used to identify a second 
gray level threshold for the initial lung segmentation. The gray level threshold is applied to 
the segmented thorax region so that the darkest pixels in the segmented thorax region remain 

1 0 "on" in the binary image that is created. 

In step 107 the anterior junction is identified if the segmented lungs are fused 
together. The presence of a single large "on" region in the binary image is used as a flag to 

- indicate that the two lungs regions are "fused" at the anterior junction. With the presence of 

an anterior junction line so established, the anterior junction line is delineated and pixels 

1 5 surrounding it are turned "off in the binary image to separate into two distinct lung regions 
what had been erroneously identified by initial gray level thresholding as a single segmented 
lung region. In step 109 the diaphragm is identified. The geometric properties of "holes" 

- 1 within the binary image caused by regions of "off pixels completely contained within larger 

"on" regions are analyzed to identify holes caused by the diaphragm. Pixels within such 

20 holes are specifically excluded from the segmented lung regions. A contour detection 
algorithm is used to identify the outer margins of the largest "on" regions in the binary 
image, and the set of all image pixels that lie within these contours (excluding pixels 
identified as diaphragm) is considered the segmented lung regions. In step 1 1 1 the 
segmented lung regions are modified by a rolling ball technique that effectively rolls a series 

25 of ball filters along the exterior aspect of the lung segmentation contours to incorporate 

pixels that may have been erroneously excluded due to initial gray level thresholding. In step 
1 13 the diaphragm is identified a second time in a second diaphragm analysis to prevent the 
rolling ball technique from incorrectly including pixels that belong to the diaphragm. 

As visualized on a CT section image, the lung regions are represented by dark (i.e., 

30 low attenuation or low CT number) regions completely surrounded by a bright (i.e., high 
attenuation or high CT number) region, which, in turn, is completely surrounded by a dark 
region (the air outside the patient). Lung segmentation proceeds by first segmenting the 
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thorax (i.e., the outer margin of the patent's body) in step 101 to eliminate from further 
consideration pixels representing air outside the patient. To achieve thorax segmentation, a 
cumulative gray level profile is constructed from the values of pixels that lie along the 
diagonal that extends from a corner of the image to the image center, as shown in Figure 2A. 
5 The shape of this profile is analyzed to identify a single gray level as a threshold value [3], 
designated by an arrow in the gray level histogram of Figure 2B. A binary image is created 
by thresholding the section image such that a pixel is turned "on" in the binary image if the 
value of the corresponding pixel in the section image has a value greater than the gray level 
threshold; all other pixels remain "off in the binary image. An eight-point contour detection 

10 scheme [27] is used to construct a contour surrounding the outermost boundary of the largest 
"on" region in the binary image (i.e., the thorax). The set of pixels in the section image that 
lie within this contour defines the segmented thorax region and is used to create a thorax 
segmentation image such that pixels within the segmented thorax region maintain their 
original value, while pixels not included within the segmented thorax region are assigned a 

15 value of 0. 

The segmented thorax region defined in this manner tends to include portions of the 
table on which the patient lies during the CT examination. The arrow in Figure 3 A points to 
a portion of the table in the image. To eliminate these pixels that represent structures 
external to the patient, each column in the thorax segmentation image is analyzed beginning 

20 at the bottom of the image (i.e., the posterior aspect). Pixels in a particular column are 

scanned until the first non-zero pixel is encountered (i.e., the first pixel within the segmented 
lung region). Subsequent pixels are examined to identify a reduction in gray level followed 
by an increase in gray level. Such a trend is assumed to represent air between portion of the 
table and the patient's body. The pixel associated with the point of maximum contrast as 

25 gray level subsequently increases is identified as the posterior margin point. The set of 
posterior margin points obtained for all image columns is smoothed to form a continuous 
margin line, and pixels that lie posterior to this margin line are eliminated from the 
segmented thorax region (i.e., assigned a gray level value of 0), as shown in Figure 3B. 

Initial lung segmentation based on gray level thresholding tends to include the trachea 

30 and main bronchi within the segmented lung regions. To ensure that these structures do not 
contribute to the segmented lung regions during the initial lung segmentation step, the 
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trachea and main bronchi are segmented and eliminated from the segmented thorax region in 
step 103. 

Figures 4A and 4B are flowcharts for explaining how the trachea and main stem 
bronchi are eliminated from the images. Referring to Figure 4A, in step 401 a seed point for 
5 trachea segmentation is automatically identified in the central region in the first (i.e., 

superiormost) thorax segmentation image. This seed point is the pixel with the lowest gray 
level (i.e., the lowest density pixel (LDP)) in a 60 by 60 pixel region centered over the center 
of mass of the thorax. In step 403 a region-growing technique [27] is performed to expand 
the identified trachea region about the seed point and thereby define the trachea region. As 

10 the gray level threshold is incremented by 5 during region growing, more pixels surrounding 
the seed point within the trachea are identified. A stopping criterion is established to halt the 
region growing process when the trachea has been adequately segmented. This stopping 
criterion is satisfied, for example, when the area of the trachea at the i* iteration is less than 5 
pixels greater than the corresponding area at the i-l th iteration. 

1 5 In step 405 the center of mass location, Cx, Cy, within the segmented trachea is 

determined. If the current image is the superiormost image, the process skips step 409 and 
returns to step 401 because step 409 compares information of two images. When the process 
returns to step 401, the next seed point is identified in the subsequent thorax segmentation 
image (i.e., the section directly below the current section). The seed point in the subsequent 

20 image is the LDP in a 15 by 20 pixel region centered over the center of mass of the trachea 
region in the previous image determined in step 405. Note that the region is larger in the y 
(vertical) direction than the x (horizontal) direction because the trachea is more likely to 
deviate in the y direction. 

Then, in step 403 region growing is again performed to segment the trachea. Next, 

25 the center of mass, Cx, Cy, of the segmented trachea is determined in step 405. 

In step 409 three conditions are checked. If any of the three conditions are satisfied, 
then the carina is assumed to have been reached in the current image, and bronchi 
segmentation begins in steps 411 and 419 (Figure 4B). If none of the three conditions in step 
409 are satisfied, then the process returns to step 401 and the seed point in the subsequent 

30 image is identified. 

The first condition is checked in step 409 by comparing (A) the horizontal location of 
the center of mass, Cx, in the current image to (B) (1) the minimum extent of the segmented 
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trachea in the x direction, minx, in the previous image and (2) the maximum extent of the 
segmented trachea in the x direction, maxx, in the previous image. If Cx < minx or if Cx > 
maxx, then the carina is assumed to have been reached in the current image, and bronchi 
segmentation begins in the current image in steps 411 and 419 (Figure 4B). 
5 The second condition is checked in step 409 by comparing (A) the gray level of the 

pixel in the current image corresponding to the center of mass of the segmented trachea in the 
previous image to (B) the region growing threshold used in step 403 for the previous image. 
If the gray level of the pixel in the current image corresponding to the center of mass in the 
previous image is more than 20 gray levels higher than the region growing threshold used for 

10 the previous image, then the carina is assumed to have been reached in the current image, and 

~ bronchi segmentation begins in the current image in steps 41 1 and 419 (Figure 4B). 

The third condition is checked by comparing the area of the segmented trachea in the 

^- current section with the area of the segmented trachea in the previous section. If the area of 

° the trachea in the current image is less than 80% of the area of the trachea in the previous 

=15 image, then the carina is assumed to have been reached in the current image, and bronchi 
segmentation begins in the current image in steps 411 and 419 (Figure 4B). 

Once the carina is determined to have been reached, two seed points (corresponding 

= : to the main stem bronchi) in each image are identified, and two region growing processes are 

performed to segment the two main bronchi in the current and subsequent images. Bronchi 

20 segmentation is not performed in images inferior to the image in which region growing first 
expands the "bronchi" into the lung parenchyma. 

Referring to Figure 4B, in step 41 1 the LDP determined in step 401 for the current 
image is designated as the seed point within the first main stem bronchus. Then, in step 413 
region growing is performed to define the first bronchus. Region growing may be performed 

25 in the same manner as in step 403. Then, in step 415 the center of mass, Clx, Cly, is 

determined for the first bronchus. If the first (i.e., superiormost) bronchi image is being 
processed, then the process returns to step 41 1 and a seed point is identified for the first main 
stem bronchus in the subsequent bronchi image. The seed points in subsequent images are 
identified as the LDP in a 30 by 20 pixel search region centered over the center of mass of 

30 the segmented first main stem bronchus determined for the previous bronchi image 
determined in step 415. The LDP search region is wider in the x direction than the y 
direction because the main stem bronchi are more likely to deviate in the x direction than the 
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y direction. Then, in step 413 region growing is performed for the new seed point 
determined in step 41 1, and in step 415 the center of mass, Clx, Cly, of the segmented 
bronchus is determined. 

The initial seed point for the second main stem bronchus is determined by finding the 
5 LDP within a 30 by 20 pixel region centered equidistant from the carina as the initial seed 
point for the first main stem bronchus. Subsequent region growing and center of mass 
determination for the second main stem bronchus in steps 421 and 423 is preferably carried 
out in the same manner as steps 413 and 415, respectively. Upon completion of step 423 for 
the second bronchus in the superiormost bronchi image, the process returns to step 419, and 
10 the seed point for the second bronchus is identified in the subsequent bronchi image. The 
i seed point is identified as the LDP in a 30 by 20 pixel search region centered over the center 

of mass, C2x, C2y, of the second bronchus in the previous image determined in step 423. 
• r Once initial segmentation is performed for the first and second bronchi in the 

: . superiormost bronchi image, two conditions are checked in step 427 for each of the first and 

15 second bronchi, after the center of mass is computed in steps 415 and 423, respectively. If 
r - ; either of the conditions is satisfied, the segmentation of the corresponding bronchus is 

terminated in step 429. If neither of the conditions is satisfied, the process returns to step 411 
J or 419. 

Referring to the first bronchus, in step 427 the first condition is checked by 
20 comparing (A) the vertical location of the center of mass, Cyl, in the current image to (B) (1) 
the minimum extent of the segmented first bronchus in the y direction, minyl, in the previous 
image and (2) the maximum extent of the segmented first bronchus in the y direction, maxyl, 
in the previous image. If Cyl < minyl or if Cyl > maxyl, then the segmentation of the first 
bronchus is terminated. The same process is applied for the second bronchus, using vales of 
25 Cy2, miny2, and maxy2, which are analogous to Cyl, minyl, and miny2, respectively. 

The second condition is satisfied if, during region growing for either the first or 
second bronchus, the number of pixels along the segmentation contour that reach the edge of 
the corresponding LDP search region is greater or equal to 20. If the second condition is 
satisfied for either bronchus, then segmentation is terminated for the corresponding bronchus. 
30 Figure 5 A is an image of the segmented trachea in a thoracic section, resulting from 

steps 403, 405, and 407. Figure 5B is an image of the segmented bronchi in another section, 
resulting from steps 413, 415, 417, 421, 423, and 425. When applied to the above-noted 17- 
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case database, this technique successfully eliminated the trachea and main bronchi from 98% 
of the segmented lung regions in which they would have otherwise been included. 

Initial lung segmentation begins for a particular section by constructing a gray level 
histogram from the pixels that lie within the segmented thorax region [3,28]. The 
5 distribution of pixels in this typically bimodal histogram is used to identify a single gray 
level as a threshold value within the broad minimum in the histogram [3]. Figure 6A is an 
exemplary gray level histogram for identifying a single gray level as a threshold value. The 
arrow in Figure 6A identifies the gray level value within the broad minimum in the 
histogram. As seen in Figure 6B, a binary image is created by thresholding the thorax 
10 segmentation image such that a pixel is turned "on" in the binary image if the value of the 
j= corresponding pixel in the thorax segmentation image has a value less than the gray level 

-v. threshold, while all other pixels remain "off in the binary image. 

» • The presence of a single "on" region that spans both sides of the resulting binary 

= i image indicates that gray level thresholding has "fused" the two lungs and that an anterior 

"15 junction line is present in the section image. Distinction between left and right lungs is often 
required for segmentation results to be useful as preprocessing for more detailed image 
analyses. According to the present invention, the single lung region separated into two 
.- 5 regions by eliminating pixels along the anterior junction line. 

Figure 7 is a flowchart of a process for separating a single lung region into two 
20 regions. In step 701, the presence of a single, large area lung, spanning both sides of the 

image is detected. The presence of a single lung are may be identified by detecting a single 
region of "on" pixels and/or by determining if a region of "on" pixels in a binary lung 
segmentation image is greater than 30% of the total area of the thorax. In step 703 a "cleft 
point" is identified. The cleft point is identified by searching the binary image for the most 
25 anterior point along the cardiac aspect of the lung regions. Figure 8A is an image of a single 
lung region with an arrow identifying the detected cleft point. Then, in step 705 the average 
pixel values along rays extending through the lung region from the cleft point are determined. 
All the rays within +/- 50 degrees of a line extending vertically from the cleft point to the 
upper edge of the lung region are analyzed for average pixel value. In step 707 the ray with 
30 the greatest average pixel value is identified as the initial anterior junction line. 

Next, in step 705, beginning with the cleft point, the set of pixels representing the 
local maximum in each row extending from the cleft point toward the anterior aspect of the 
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lung region is determined. This search is performed within +/- 10 pixels of the initial 
anterior junction line along each row of the initial anterior junction line. These local 
maximum pixels are designated the anterior junction line in step 711. Then, in step 713 the 
local maximum pixels, along with one pixel on either side of each local maximum pixel, are 
5 turned "off in the binary image. As a result, two distinct regions within the binary image are 
created, as shown in Figure 8B. When applied to the 17-case database, this technique 
accurately delineated the anterior junction line in all section images in which an anterior 
junction line was present (100% accuracy). 

An eight-point contour detection scheme [27] is used to construct contours 
10 surrounding the outermost boundaries of the two largest "on" regions (i.e., the lungs) in the 

!=, binary image in Figure 6B. The sets of pixels in the section image that lie within these 

contours define the segmented lung regions and are used to create a lung segmentation image 

"~ such that pixels within the segmented lung regions maintain their original value, while pixels 

not included within the segmented lung regions are assigned a value of 0 (Figure 1 1 A). 

? 1 5 As evident from Figure 6B, the binary images that result from gray level thresholding 

tend to contain "holes" of "off pixels that are completely surrounded by "on" pixels. These 
holes result from denser (i.e., brighter) structures contained within the lung regions that have 
gray levels greater than the gray level threshold for initial lung segmentation. Consequently, 
the pixels corresponding to these denser structures remain "off in the binary image. In most 

20 instances, these structures represent vessels that are part of the anatomic lungs. Since the 
contouring scheme considers a segmented lung region as all pixels within the outermost 
boundary of an "on" region in the binary image, these dense vessels are correctly included 
within the segmented lung regions. However, the diaphragm, which is not part of the lung, 
often results in a similar hole in the binary image (e.g., as shown by the arrow in Figure 

25 12 A). 

According to the present invention, each hole in a binary image is identified, and 
features such as area and circularity are computed in order to prevent the inclusion of pixels 
that belong to the diaphragm. In this manner, holes caused by the diaphragm may be 
identified, and the corresponding pixels may be excluded from the segmented lung regions 
30 using the processes described in Figures 13A and 13B, for example. 

The segmented lung regions based on gray level thresholding tend to exclude dense 
structures along the edges of the lung regions such as juxta-pleural nodules and hilar vessels. 
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Figure 9 is an image of the initial lung segmentation contours with arrows identifying a juxta- 
pleural nodule and hilar vessels. 

A rolling ball algorithm is applied to properly include within the segmented lung 
regions dense structures along the edges of the lung regions [14,29]. A circular filter (the 
5 "ball") is constructed and is "rolled" along the lung segmentation contours by successively 
identifying that pixel along the ball's circumference with a tangential slope that matches the 
slope of the current contour point. The filter is then positioned to align the selected ball 
circumference pixel with the contour pixel. If an indentation of the proper scale is 
encountered, the ball will overlap the contour at some contour point other than the point of 
10 contact used to place the filter. This overlap point and the point of contact define endpoints 
1 ! of the indentation. Linear interpolation is then used to create new contour points that connect 

I these endpoints, effectively bridging the gap in the contour and eliminating the indentation. 

Figure 10A is an image of two lung segmentation contours. The boxed region in 
Figure 10A is shown expanded in Figure 10B with a rolling ball filter overlapping the point 
" = 15 of contact used to place the filter and one other point (i.e., the two points of contact). The 
result of linear interpolation between the two points of contact is shown as a dark line 
connecting the two points. The newly encompassed image pixels are added to the segmented 
lung regions. 

Figure 1 1 A is an image of the lung contours before application of the rolling ball 
20 filter, and Figure 1 IB is an image of the lung contours after application of the rolling ball 

filter. An iterative, multi-scale approach is used in which balls of different radii are applied 
in succession to rectify indentations of different dimensions. The technique can also be 
applied in three-dimensions by rolling a spherical filter along the external aspect of the set of 
segmented lung regions from all sections (surfaces) in a case considered as a complete 
25 volume. 

One potential pitfall of the rolling ball algorithm is that it may force inclusion of the 
diaphragm, as shown in Figure 12B, even when the diaphragm has been properly excluded by 
gray level thresholding. To rectify this problem, the present invention computes geometric 
features of indentations identified by the rolling ball algorithm. Based on the values of these 
30 features, the rolling ball algorithm is prevented from bridging indentations caused by the 
diaphragm. 
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Figures 13 A and 13B are flowcharts of processes for identifying pixels that belong to 
the diaphragm and excluding such pixels from the segmented lung regions. The process of 
Figure 13A is applied in step 109 of Figure 1. In step 1301, all holes that exist within the 
binary image created in step 105 are identified and labeled. In step 1303 at least one 
5 geometric feature of each hole (e.g., the area of each hole) is computed. Then, in step 1305 it 
is determined whether the geometric features of each hole exceed predetermined thresholds. 
Those holes having geometric features that exceed the predetermined values are identified as 
corresponding to the diaphragm. For example, if the area of the hole is greater than a 
predetermined threshold value of 706 mm 2 , then the hole is identified as part of the 

10 diaphragm. Preferably, geometric feature processing is performed only on the bottom half of 
the images, since the diaphragm is in the lower part of the thorax. In step 1307 the pixels 

_ forming holes identified as corresponding to the diaphragm are excluded from the segmented 

- lung region. 

The process of Figure 13B is applied during the application of the rolling ball 

"1 5 algorithm in step 1 1 3 of Figure 1 . In step 1 309 the rolling ball algorithm is performed to 

identify indentations along the periphery of the segmented lung region, as described above. 
In step 1311 at least one geometric feature of each indentation is determined. In step 1313 
the geometric features of each indentation are compared against predetermined values to 
determine whether the geometric features of each indentation exceed predetermined values. 

20 Those indentations having geometric features that exceed the predetermined values are 

identified as corresponding to the diaphragm. In a preferred embodiment, the predetermined 
threshold value for two geometric features must be exceeded in order for an indentation to be 
identified as part of the diaphragm. The first feature is the number of pixels in the 
indentation (measured between the connecting points of the rolling ball) divided by the total 

25 number of pixels forming the entire lung contour of the corresponding lung. The threshold of 
the first feature is 0.20. The second feature is the compactness of the indentation. 
Compactness is equal to (A) the area encompassed by the indentation and the line segment 
connecting the two connecting points of the rolling ball divided by (B) the area of a circle 
having a circumference equal to the perimeter of the indentation. The perimeter of the 

30 indentation is the length of the line segment connecting the two connecting points of the 
rolling ball plus the length of the contour between the two connecting points. 
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Then, in step 13 15 the rolling ball algorithm is prevented from bridging indentations 
corresponding to the diaphragm in order to prevent the rolling ball algorithm from including 
indentations corresponding to the diaphragm within the segmented lung regions. 

Figure 14 is a block diagram of a system for segmenting lung regions in thoracic CT 
5 images. The blocks in Figure 14 correspond to program modules, circuits, and/or 

mechanisms configured to implement the method(s) described above. CT scans of an object 
are obtained from an image acquisition device 1401 and input to the system. Each image is 
stored in memory 1403. The image data of each section image from a particular CT scan is 
first passed through the cumulative gray level profile circuit 1405 and then to the gray level 
10 profile analysis circuit 1407 for gray level threshold selection. The image data along with the 
I) gray level threshold value are passed through the gray level thresholding circuit 1409 and 

modified by passing through the table detection circuit 1411. The data are then passed 
j ; through the contour construction circuit 1413. The image data are passed through the trachea 
j: and main stem bronchi detection circuit 1415 prior to being sent through the gray level 
i5 histogram circuit 1417. The output from the gray level histogram circuit are sent to the 

histogram analysis circuit 1419 for gray level threshold value identification. The image data 
along with the gray level threshold value are passed through the gray level thresholding 
=. circuit 1409 and modified by anterior junction circuit 1421. Next, the image data is passed 

through the first diaphragm detection circuit 1423 and then through the contour construction 
20 circuit 1425. The contours are modified through the rolling ball circuit 1427 which includes 
the second diaphragm detection circuit 1429. In the superimposing circuit 1431 the results 
are either superimposed onto images, stored in file format, and/or output in text format. The 
results are then displayed on the display system 1433 after passing through a digital-to- 
analog converter 1435. 

25 This invention conveniently may be implemented using a conventional general 

purpose computer or microprocessor programmed according to the teachings of the present 
invention, as will be apparent to those skilled in the computer art. Appropriate software can 
readily be prepared by programmers of ordinary skill based on the teachings of the present 
disclosure, as will be apparent to those skilled in the software art. The invention may also be 

30 implemented by the preparation of application specific integrated circuits or by 

interconnecting an appropriate network of conventional component circuits, as will be readily 
apparent to those skilled in the art. 
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Figure 15 is a schematic illustration of a computer system for segmenting lung 
regions in CT scans. A computer 1500 implements the method of the present invention, 
wherein the computer housing 1502 houses a CPU 1506, memory 1508 (e.g., DRAM, ROM, 
EPROM, EEPROM, SRAM, SDRAM, and Flash RAM), and other optional special purpose 
5 logic devices (e.g., ASICs) or configurable logic devices (e.g., GAL and reprogrammable 
FPGA). The computer 1500 also includes plural input devices, (e.g., a keyboard 1522 and 
mouse 1524), and a display card 1510 for controlling monitor 1520. In addition, the 
computer 1500 further includes a floppy disk drive 1514; other removable media devices 
(e.g., compact disc, tape, and removable magneto-optical media); and a hard disk 1512, or 

10 other fixed, high density media drives, connected using an appropriate device bus (e.g., a 
SCSI bus, an Enhanced IDE bus, or a Ultra DMA bus). Also connected to the same device 
bus or another device bus, the computer 1500 may additionally include a compact disc 
reader, a compact disc reader/writer unit or a compact disc jukebox. 

As stated above, the system includes at least one computer readable medium. 

1 5 Examples of computer readable media are compact discs, hard disks, floppy disks, tape, 
" ' magneto-optical disks, PROMs (EPROM, EEPROM, Flash EPROM), DRAM, SRAM, 

SDRAM, etc. Stored on any one or on a combination of computer readable media, the 
present invention includes software for controlling both the hardware of the computer 1500 
and for enabling the computer 1500 to interact with a human user. Such software may 

20 include, but is not limited to, device drivers, operating systems and user applications, such as 
development tools. Such computer readable media further includes the computer program 
product of the present invention for performing the inventive method described above. The 
computer code devices of the present invention can be any interpreted or executable code 
mechanism, including but not limited to scripts, interpreters, dynamic link libraries, Java 

25 classes, and complete executable programs. Moreover, parts of the processing of the present 
invention may be distributed for better performance, reliability, and/or cost. For example, an 
outline or image may be selected on a first computer and sent to a second computer for 
remote diagnosis. 

Obviously, numerous modifications and variations of the present invention are 
30 possible in light of the above teachings. For example, the outline of the nodules may be 

extracted using any available automated technique, rather than manually. It is therefore to be 
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understood that within the scope of the appended claims, the invention may be practiced 
otherwise than as specifically described herein. 
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